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We show that time dependent gravitational potential can be directly detected from the cosmic 
microwave background (CMB) anisotropies. The signature can be measured by cross-correlating 
the CMB with the projected density field reconstructed from the weak lensing distortions of the 
CMB itself. The cross-correlation gives a signal whenever there is a time dependent gravitational 
potential. This method traces dark matter directly and has a well defined redshift distribution of 
the window projecting over the density perturbations, thereby avoiding the problems plaguing other 
proposed cross-correlations. We show that both MAP and Planck will be able to probe this effect for 
observationally relevant curvature and cosmological constant models, which will provide additional 


constraints on the cosmological parameters. 

PACS numbers: 98.80.Es,95.85.Nv,98.35.Ce,98.70 

It is widely accepted that cosmic microwave back¬ 
ground (CMB) anisotropies offer a unique environment to 
study cosmological models. The anisotropies were gen¬ 
erated predominantly during recombination at redshift 
z ~ 1100, when the universe was still in a linear regime 
and the physics at eV energy scale was simple. This 
allows one to make robust predictions for various cos¬ 
mological models, which can be compared to an increas¬ 
ing number of observations. However, some degeneracies 
between cosmological parameters remain even for future 
satellite missions and these are being further expanded as 
new parameters are being introduced. The degeneracies 
are particularly severe between various components af¬ 
fecting the expansion of the universe, such as curvature, 
cosmological constant or any other term with a more gen¬ 
eral equation of state |jj. Other cosmological tests must 
therefore be used to break these degeneracies. 

It has been pointed out that in a universe where matter 
density does not equal critical density the gravitational 
potential is changing with time, which produces a signif¬ 
icant component to the CMB on large scales j|. This 
effect is generated at late times and since the gravita¬ 
tional potential is related to the density field through 
the Poisson’s equation, the effect can also be looked for 
by cross-correlating CMB with another tracer of density 
field H. Unfortunately, no clean density map out to high 
redshift exists on large scales. The X-ray background has 
been suggested as a possible tracer of large scale structure 
out to z ~ 3 — 4, but the uncertainties associated with 
the redshift distribution of the sources, relation between 
X-ray light and underlying mass and the Poisson fluctua¬ 
tions from the nearby sources make this test inconclusive 



Recently we developed a method to reconstruct the 
projected density field out to recombination directly from 
the CMB anisotropies |]]. The method is based on the 


gravitational lensing effect, which distorts the pattern of 
CMB anisotropies j?J. Although the signal to noise for 
individual structures from such a reconstruction is small, 
averaging over independent patches of CMB reduces the 
noise and extracts the signal in a statistical sense. We 
were able to show that this allows one to extract the 
power spectrum of density perturbations with high accu¬ 
racy over two decades in angular scale [||. 

In the present paper we use the reconstructed pro¬ 
jected density field to cross-correlate it with the CMB 
itself. If there is a component to the CMB from the time 
evolving gravitational potential then it should correlate 
with the projected density field. Most of the signature 
comes from large angular scales, so we first generalize 
the method developed in Q to all sky. Because the 
small scale CMB anisotropies were generated uniquely 
during recombination, the weighting of density pertur¬ 
bations as a function of redshift in the projection is well 
defined. Moreover, gravitational lensing effect depends 
on the dark matter distribution in the universe, so no 
assumption of how light traces mass is necessary. This 
avoids the shortcomings of cross-correlations with X-ray 
and other tracers of large scale structure mentioned in 
[|J. In addition, the projected density field is sensitive 
to matter distribution out to a very high redshift and 
allows one to test the models where the time dependent 
potential is generating anisotropies at higher redshifts, 
such as the curvature dominated models |yj. Here we 
compute the expected signal to noise of future CMB mis¬ 
sions for cosmological constant and curvature dominated 
models, using both MAP and Planck satellite character¬ 
istics. Although we limit to these two families of models 
we note that other models, such as those with more gen¬ 
eral equation of state, would also produce a signature 
that one could look for. 

To reconstruct the projected density field we consider 
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the symmetric tensor of products of temperature deriva¬ 
tives transverse to direction n |(J 

n ab = —T :a T, b ( 1 ) 

0\s 

where T :a , T :b are covariant derivatives of T with respect 
to the coordinate basis in the tangent space of direction 
n, here defined with polar coordinates (8, (f>) and g ab is 
the metric on the sphere. We defined as so that in the 
absence of lensing the average over CMB (H ab )cMB = 
\g ab . The tensor can be decomposed into the trace and 
traceless component as H ab = ^(1 — S)g ab — 'Hah- From 
the traceless tensor H ab one may define two rotationally 
invariant quantities 

£ = l^~ 2 Ha b b B = \v~ 2 H a a b c e a c , (2) 

where e“ is the completely antisymmetric (Levi-Civitta) 
tensor and V~ 2 is the inverse Laplacian on the sphere. 

In the presence of lensing the average of H ab becomes 

§ 

C Hab)CMB = ^2 _ K ^j 9a b — lab (3) 

where k and ^ ab are the convergence and the shear com¬ 
ponents of the symmetric shear tensor & ab , defined as 
the covariant derivative of the displacement field on the 
sphere ]fj], which encodes the information on the grav¬ 
itational lensing effect. All rotationally invariant quan¬ 
tities can be decomposed on a sphere into spherical har¬ 
monics, X(n) = im. a x,imYim(n), where X stands for 
T, k, S, £ or B. From equation |Q we find (as,i m } = 
Msio-K,im, with Msi = 2. The multipole moments of 
the scalar field £ average to ( a,£,i m ) = M£ia K j m , where 
Msi = 2(7 + 2 )(l — 1 )/l(l + 1) Q, while the average 
of the pseudo scalar field B identically vanishes in the 
large scale limit, ( as,im) = 0, because gravitational po¬ 
tential from which shear is generated is invariant under 
the parity transformation. Convergence k can thus be 
reconstructed in two independent ways from S and £, 
while the third quantity B serves as a check for possible 
systematics. Note that since convergence is expressed 
as a quadratic quantity of T its cross-correlation with T 
gives a non-vanishing 3rd moment. This means it can 
also be looked for using bispectrum, which is a method 
independently proposed by (To| . 

Convergence can be written as a projection of grav¬ 
itational potential QQ k = fg° g(x, Xo)V 2 (f>(x)dx- 
Here xo is the comoving radial coordinate at recombi¬ 
nation and g is the radial window, defined as g{XiX o) = 

■ It is a bell shaped curve symmetric around 
x/2 and vanishing at 0 and xo- Here r(x) is the comoving 
angular diameter distance, defined as AT -1 / 2 sinAT 1/,2 x, 
X, (—A') -1 / 2 sinh(—AT) 1 / 2 x for K > 0, K = 0, K < 0, 


respectively, where K is the curvature. Curvature can 
be expressed using the present density parameter fio and 
the present Hubble parameter Hg as K = (Hq — l)Hg . In 
general Ho consists both of matter contribution fi m and 
cosmological constant term Ha- 

The angular power spectra are defined as C xx = 
aY+rYjm a *x im a x' ,im- Their ensemble averages are given 

by |! 

Cf*'= (4 tt) 2 / f3 2 d(3P((3)A xl (/3, To )A X ' l ((3,To) 
A X i(f3, t 0 ) = / 0 TO dr$jg(T 0 - t)S x {/3, t), (4) 

where are the ultra-spherical Bessel functions and 

P(/3) is the primordial power spectrum. Equation |]only 
applies to flat and open universes, whereas for the closed 
universe the eigenvalues of the Laplacian are discrete so 
the integral over (3 is replaced with a sum over K~ ? (3 = 
3,4, 5.... The source for temperature anisotropies St is a 
combination of several terms. These can be decomposed 
into terms generated during recombination, which consist 
of Sachs-Wolfe term, Doppler term, intrinsic anisotropy 
term and anisotropic stress term, and late time term gen¬ 
erated by the time dependent gravitational potential (so- 
called integrated Sachs-Wolfe term or ISW). The latter is 
only important for low multipole moments. The full form 
of St can be found in 0- The source for convergence 
is S K = g(3 2 (j). It is also useful to define the correlation 
coefficient Corr^ K = Cf K /(C' J TT C' i KK ) 1 / 2 , which is the rel¬ 
evant quantity for the estimation of signal to noise. 

Using the above expressions one can compute Cf T , 
Cf K and Cj K for any cosmological model. We performed 
numerical calculations using a modified version of CMB- 
FAST code (u|. We have verified that agrees with 
previous calculations which were done in the small scale 
limit jn|, as well as with the alternative all-sky expres¬ 
sions given in |H|. The results for (Corrf K ) 2 are shown 
in figure |l] for a representative set of models. The corre¬ 
lation coefficient in a cosmological constant model is sub¬ 
stantial only for very low l and rapidly drops on smaller 
scales. In a curvature dominated model with the same 
Q rn the correlation coefficient is larger to start with and 
also drops less rapidly with l. This indicates that one 
will be able to set stronger limits on curvature than on 
cosmological constant, which is confirmed below with a 
more detailed analysis. The reason is that in a cosmolog¬ 
ical constant model the gravitational potential changes 
only at late times (z'S 1 for reasonable values of H m ), 
while in a curvature dominated model potential changes 
also at higher redshift. This leads to two effects. First is 
that in a cosmological constant model ISW is compara¬ 
ble to other terms only for the lowest multipoles, while 
in a curvature model ISW dominates up to higher l |5J. 
Second is that the window g peaks at relatively high red- 
shift x ~ 3 and so is able to pick up correlations with 
ISW from open universe better than that from cosmo¬ 
logical constant universe. We have also calculated the 
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correlation coefficient for a flat model. It is small for 
all l, (Corr^ K ) 2 < 10 —3 , demonstrating that correlations 
with fluctuations generated at recombination are negli¬ 
gible and the cross-correlation is sensitive to the time 
dependent gravitational potential only. 



FIG. 1. Square of correlation coefficient (Corr^ K ) 2 as a 
function of l is shown for open = 0.2 and = 0.4 
models (upper and lower dashed curve) and for cosmological 
constant = 0.2 and = 0.4 models (upper and lower 
solid curves). 

We now address the question of signal detectability 
with the future CMB missions. We continue to work 
in multipole moment space and assume we have all sky 
expansion, which allows us to decouple between differ¬ 
ent m and l multipole moments. The generalization to 
incomplete sky coverage can be approximated by insert¬ 
ing appropriate factors of sky coverage fraction f s j- y in 
the final expressions. Given two random fields T and 
W (where W stands for S or £) we want to develop a 
test that maximizes the signal in the presence of cor¬ 
relations against the null hypothesis that there are no 
correlations. The term that quantifies the correlations 
is the product between the two fields Xj^ = T* m Wi m 
(here and below average with the complex conjugate 
(Tj^ n VVi m +Ti m yVi rn )/2 is implied). Its expectation value 
under the null hypothesis of pure noise is (Xj ^)o = 0, 
because the function entering this expression is a three- 
point function of T, which vanishes both for intrinsic 
fluctuations and for detector noise under the Gaussian 
assumption. The alternative hypothesis is that of sig¬ 
nal which gives {X'^ l ) 1 = {T^ m Wi m ) = M W iC? K . The 
variance under the null hypothesis is 

/ V W vW' \ /yW\ / vW'\ 

X^lm^lm /0 \^-Im)0 \^-lm /O 

(Q tt + N? T )(M wl M wn Cr + (5) 

where -/V ; TT and N^ w are the noise power spectra for 
CMB anisotropies, S, £ or their cross-term, respectively. 


Both S and £ contribute information. If they 
are uncorrelated then the information contents can be 
added independently, otherwise the covariance matrix 
Cov{X'^ l X^) has to be diagonalized first. The CMB 
term Cj T + N^ T is the same for all matrix elements and 
can be computed using MAP and Planck noise charac¬ 
teristics. For these CMB missions detector noise on large 
scales will be negligible, hence A) TT <C Cj T . The domi¬ 
nant source of noise in S or £ on large scales are the CMB 
anisotropies. The noise terms involve integrals over the 
CMB power spectrum and can be computed numerically 
using the expressions given in Q . The results are shown 
in figure || both for MAP and Planck. They show that 
on large scales the CMB noise power spectrum has ap¬ 
proximately white noise shape, N^ w ~ const. At low l 
Nf £ ~ N? s /2 > Nf s , C? K . Therefore, noise dominates 
over the signal and the latter can only be extracted in 
a statistical sense by averaging over multipole moments. 
Because the off-diagonal term ( £* m Si m )o is much smaller 
than the two diagonal terms the covariance matrix is 
nearly diagonal and the information from S and £ can 
be added independently, with £ contributing twice the 
amount of information than S. Note also that Planck 
has a factor of 5 better sensitivity than MAP. 



FIG. 2. Power spectra for noise Nf s (long dashed), N ££ 
(short dashed), N £S (dotted) both for Planck (lower curves) 
and MAP (upper curves) for cosmological constant model 
with fi m = 0.3. Also shown is power spectrum of conver¬ 
gence 4 Ci* (solid) for the same model, normalized to ag = 1. 

We now want to combine the signal to noise from 
different multipole moments to maximize the overall 
signal to noise. To do this we add up the products 
weighted with some yet to be determined weights a;, 
X = ^2m i a iXim- Since the moments are uncorrelated 
the expectation value and variance are 

(A-)i =^(2Z + 1 )aiM Wtl C? K 

l 
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(x 2 ) 0 = +1 W(q tt + Nr)(M^cr+ n^ w ), 

i 

while the null hypothesis mean remains (X)o = 0. We 
want to maximize S/N = ((X)i — (X)q)/((X 2 )o) * 1 11 ^ 2 with 
respect to a;. Taking derivatives with respect to ai and 
setting the expression to 0 we find ai = C) Tk /(Cj T + 
Nf T )(XI^ vl Cf K + N^ vw ). The overall signal to noise is, 
combining the information from S and £ 

S_ 

N 


(7) 


fsky 


(21 + l)(Corr/ K ) 


l W=£,5 (.1 + JyTt j(,l 


Jk\2 




N ww ^ 
Tjir, t) 




where we inserted f s k y to account for the fact that the 
number of multipoles will be smaller if only some fraction 
of the sky will be measured. We have expressed the sig¬ 
nal in terms of number of standard deviations above the 
noise. To express it in terms of confidence limits a more 
detailed analysis with Monte Carlo simulations is needed 
H, but in the limit where many multipole moments con¬ 
tribute to the signal the usual Gaussian confidence limits 
as a function of number of standard deviations apply. 
The above expression shows that if correlation is unity 
and noise is negligible then each multipole moment con¬ 
tributes one degree of freedom and the signal to noise is 
as expected N d ' o{ , where iVdof is the number of degrees 
of freedom. Decorrelation and/or noise decrease the ef¬ 
fective number of degrees of freedom. 

Using above expressions we find S/N = 8 for £l m = 
0.4 open model and S/N = 13 for £l m = 0.2 open 
model, both for Planck noise and beam properties us¬ 
ing fsky = 0.7. Corresponding numbers for MAP are 
3.5 and 7. Both MAP and Planck will thus be able to 
usefully constrain open models with f i m < 0.4, which 
spans the range of currently favored values of fl m . For 
cosmological constant model the numbers are somewhat 
lower, Planck gives S/N = 3 and 6 for fl m = 0.4 and 
0, m = 0.2, respectively, while corresponding MAP num¬ 
bers are 1 and 2. A positive detection in these mod¬ 
els can therefore only be obtained with Planck, unless 
is very low. One can use the absence or presence of 
cross-correlation to put constraints on the models. Any 
detection of the signal with MAP will for example be 
more easily explained in terms of curvature models than 
with cosmological constant models, while absence of the 
signal in Planck will certainly rule out all curvature mod¬ 
els of interest, as well as put strong constraints on cos¬ 
mological constant models. Within the context of more 
specific models, such as the family of CDM models, one 
can use the cross-correlation to break the degeneracies 
present when only the CMB power spectrum constraints 
are used. The well-known degeneracy between curva¬ 
ture and cosmological constant can for example be bro¬ 
ken using this cross-correlation. A more detailed analy- 


(6)fis which includes temperature, polarization and conver¬ 
gence information will be presented elsewhere (h]]. Note 
that the theoretical limit for signal to noise can be ob¬ 
tained by assuming k is perfectly known and is given by 
S/N = + l)(Corrf K ) 2 . This gives S/N about a 

factor of 2 higher than our results for Planck above. 

Finally, we should mention that signal should be con¬ 
sistent with the null hypothesis (pure noise) for the B 
field in the large scale limit. Any evidence against that 
would be a sign of a systematic effect present in the data. 
This test provides a useful overall check of the method. 
Another useful test would be cross-correlating £ and S 
with the polarization CMB map. Since ISW does not 
contribute to the latter the result should again be consis¬ 
tent with zero. The straightforward interpretation and 
many consistency checks make the here proposed method 
one of the most promising ways to determine cosmolog¬ 
ical parameters and should provide further incentive for 
high sensitivity all-sky CMB experiments. 
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